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Abstract 

We compute the electric-current susceptibility x °f hot quark-gluon matter in an 
external magnetic field B. The difference between the susceptibilities measured in 
the directions parallel and perpendicular to the magnetic field is ultraviolet finite 
and given by x ~ X ± = VTN C £\ 9/|9/-B|/(2vr 2 ), where V denotes the volume, T 
the temperature, N c the number of colors, and qf the charge of a quark of flavor /. 
This non-zero susceptibility difference acts as a background to the Chiral Magnetic 
Effect, i.e. the generation of electric current along the direction of magnetic field in 
the presence of topological charge. We propose a description of the Chiral Magnetic 
Effect that takes into account the fluctuations of electric current quantified by the 
susceptibility. We find that our results are in agreement with recent lattice QCD 
calculations. Our approach can be used to model the azimuthal dependence of charge 
correlations observed in heavy ion collisions. 



1 Introduction 



The strong CV problem - the absence of V and CP violation in strong in- 
teractions - still stands as one of the fundamental puzzles of contemporary 
physics. The puzzle stems from the existence of vacuum topological solutions 
in Quantum Chromo-Dynamics (QCD) pQ. Due to the axial anomaly [2] these 
topological solutions induce the non-conservation of flavor-singlet axial cur- 
rent. The resulting time-dependence of the axial charge leads to the picture 
of the physical QCD vacuum representing a Bloch-type superposition of an 
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infinite number of topologically distinct but degenerate in energy sectors con- 
nected by the tunneling transitions - so called vacuum" [3]. Instead of 
keeping track of this vacuum structure explicitly one can instead equivalently 
reproduce its effect by adding to the QCD Lagrangian a new "#-term" , 9 ■ Q, 
where Q = (g 2 /16ir 2 ) tvF^F^ is the density of topological charge. Unless 
9 = 0, this term explicitly breaks V and CV symmetries of QCD. 

Among the effects induced by the #-term is a generation of the (V- and CP- 
odd) electric dipole moments (e.d.m.'s) of hadrons. The current experimental 
upper bound on the neutron's e.d.m. is \d n \ < 2.9 x 10 -26 e • cm Since 
inducing a non-zero electric dipole moment requires flipping the chirality of 
the quark that is achieved by the quark mass m q insertion, on dimensional 
grounds one expects \d n \ ~ (em q /m%) 9 ~ 1CT 16 9 e-cm where rriN is a typical 
hadronic scale that we have chosen as a nucleon mass. This rough estimate 
combined with the experimental measurement of neutron's e.d.m. leads to 
the bound \9\ < lCT 10 . (A more careful analysis yields a somewhat tighter 
bound, \9\ < 0.7 x 10~ n [5].) The mechanism responsible for the unnatural 
smallness of the parameter 9 has not been established yet - this represents the 
strong CV problem. One appealing explanation (that however has not been 
confirmed yet) promotes 9 into a dynamical axion field [5117] emerging as a 
Nambu-Goldstone boson of an additional chiral symmetry [8]; for a review, 
see Ref. [9]. 

The strength of topological charge fluctuations in QCD vacuum is quantified 
by topological susceptibility that is defined as a second derivative of the QCD 
partition function with respect to 9. On the other hand, the dependence of 
the partition function on 9 at low energies is governed by chiral symmetry; 
as a result, the topological susceptibility can be expressed in terms of and 
(qq) [TOlfTT] . The phase structure associated with finite 9 (especially at 9 « 
7r) [12] can also be investigated in the framework of chiral effective models |13j . 
The fluctuations of topological charge affect the mass spectrum and other 
properties of hadrons. 

At high temperature and at weak coupling, topological fluctuations in QCD 
matter are enhanced due to the real-time "sphalerons" (akin to the thermal 
activation processes) [THI15|I16|I17II18] . At temperatures not far from the decon- 
finement transition temperature T c the physics is essentially non-perturbative, 
and the description of real-time dynamics has to rely on numerical simulations 
or models. It was envisioned early on that V- and CP-odd condensates [19] 
may develop locally in dense hadronic matter [20J. It is thus conceivable that 
V and CV symmetry may be violated locally [21J. 



The weak interactions induce CV violation through the phases in the CKM 
matrix, but the resulting neutron's electric dipole moment is very small, |<i n | ~ 
10~ 32 e • cm. 
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An explicit metastable solution describing a P- and CP-odd "bubble" in hot 
QCD matter at T ~ T c has been found by using a chiral lagrangian descrip- 
tion [22J. This metastable solution describes a local domain filled with the rj' 
condensate (or, equivalently, characterized by a locally non- vanishing 9); it is 
somewhat analogous to the disoriented chiral condensate [23J that is however 
unstable even classically. It has been proposed [24J that in heavy ion collisions 
P-odd bubbles would induce certain P-odd correlations in pion momenta, but 
the experimental study of these correlations appeared challenging experimen- 
tally [25125] . 

However some time ago it was proposed that the presence of magnetic field 
and/or angular momentum in heavy ion collisions opens new possibilities for 
the observation of P- and CP-odd effects. Specifically, it was found that in 
the presence of magnetic field and/or angular momentum the fluctuations of 
topological charge can be observed directly since they lead to the separation of 
electric charge along the axis of magnetic field due to the spatial variation of 
the topological charge distribution [271l28f29] and to the generation of electric 
current due to the time dependence of the topological charge density 2?) 30113 lj. 
i.e. the "Chiral Magnetic Effect" (CME). 

The experimental observable that is sensitive to this locally P- and CP-odd 
charge separation in heavy ion collisions has been proposed in Ref. [32]. The 
preliminary data from STAR Collaboration at Relativistic Heavy Ion Col- 
lider at Brookhaven [3"3"f3"l] indicated the presence of the charge-dependent 
azimuthal correlations, with the magnitude consistent with the early rough 
estimate [27]. Recently, STAR Collaboration presented the conclusive obser- 
vation of charge-dependent azimuthal correlations [35J possibly resulting from 
the (P- and CP-odd) charge separation with respect to the reaction plane. 
Because the sign of charge asymmetry is expected to fluctuate from event to 
event, these measurements are performed on the event-by-event basis. The 
experimental observable [52] measures the strength of charge asymmetry fluc- 
tuations; even though it is sensitive to parity-violating effects, it is P- and CP- 
even and so "conventional" backgrounds have to be carefully studied. None 
of the existing event generators (such as MEVSIM, UrQMD, and HIJING) 
can reproduce the observed effect [35J although the search for other possible 
explanations of course has to continue; see e.g. [3BJ. At the same time, it is of 
paramount importance to establish a firm theoretical framework allowing for 
a quantitative study of P- and CP-odd effects in heavy ion collisions. 

The STAR result has already excited significant interest, and a number of 
recent theoretical studies address in detail the physics of the Chiral Magnetic 
Effect. The first lattice study of CME has been performed by the ITEP lattice 
gauge group [37] in quenched QCD. Recently, the Connecticut group [38] 
performed the first study of CME in full QCD with (2 + 1) light-quark flavors 
in the domain-wall formulation. These lattice studies provide an important 
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confirmation of the existence of CME, but also highlight a need for a quanti- 
tative theoretical understanding of the involved non-perturbative phenomena. 
In this paper we will attempt to reproduce some of the lattice results in an 
analytical approach. 

The behavior of CME at strong coupling in the Sakai-Sugimoto model and 
related theories has been explored through the AdS/CFT correspondence in 
Refs. [39ll40|41|42|43j . Some of these results at present are under discussion 
- for example, while Yee finds in Ref. [10] that the magnitude of CME at 
strong coupling is not modified relative to the weak coupling case, the authors 
of Ref. [H] argue that the effect disappears in the strong coupling limit. In 
Ref. [H] the CME at low temperatures has been studied using the instanton 
vacuum model. The electric dipole moment of QCD vacuum in the presence 
of external magnetic field has been evaluated in Ref. [15] using the chiral per- 
turbation theory. Extensive studies of the local V- and CV violation in hot 
hadronic matter and of the influence of magnetic field on the phase transitions 
have been performed in Ref. [16]. The properties of hadronic matter in external 
magnetic fields have attracted significant interest recently [4Tll4^] . The ana- 
lytical [30] and numerical [4"5f5T)] evaluations of the strength of magnetic field 
produced in heavy ion collisions yield the values of the order of \eB\ ~ rn* for 
RHIC energies, making it possible to study the interplay of QCD and QED 
phenomena. The violation of parity in cold and dense baryonic matter has 
been investigated in Refs. |51j . A closely related problem is the ^-dependence 
in hadronic matter at finite baryon and isospin densities [52] as well as at high 
temperatures [53] . 

Let us now discuss the physics of CME in more detail. First, note that a con- 
stant (homogeneous in space and time) 6^0 cannot induce the CME. This 
is because the #-term in the QCD lagrangian represents a full divergence and 
thus cannot affect the equations of motion. The CME electric current can be 
induced only when 9 changes in time from a finite value in the metastable state 
toward zero in the ground state; a spatially inhomogeneous 9 distribution in- 
duces the electric dipole moment [28.29J. To give a quantitative description of 
the CME, in Ref. [31] we used the chiral chemical potential /X5 that is propor- 
tional to the time derivative of 9: ^5 = 8q9/ (2Nf). The anomaly relation allows 
us to find an exact expression for the induced current which is proportional 
to B and ^5. In analogy to the ordinary relation between the current and the 
electric field through the electric conductivity, the coefficient of the current in 
response to B can be called the "chiral magnetic conductivity." In Ref. [51] 
the chiral magnetic conductivity was computed as a function of the energy 
and the momentum at weak coupling. The chiral magnetic conductivity was 
then evaluated by holographic methods also at strong coupling [40J. 

The experimental observable studied by STAR [35|55] measures the strength of 
event-by-event fluctuations of charge asymmetry relative to the reaction plane, 
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i.e. along the direction of magnetic field (and of orbital momentum). Therefore 
the quantity that needs to be evaluated theoretically is the correlation function 
of electric charge asymmetry; under some reasonable assumptions it can be 
related to the susceptibility of the CME electric current. This quantity has 
also been computed recently in the lattice QCD simulation [37f38] . 

In this paper we will compute the electric-current susceptibility in an analyt- 
ical approach. We hope that this calculation will serve as a step in a quanti- 
tative theoretical understanding of both experimental and lattice results. In 
our study we consider explicitly only the quark sector; we justify this by the 
absence of perturbative corrections to the axial anomaly that is driving the 
CME. This is justified at weak coupling - in other words, the CME current is 
not perturbatively renormalized. Whether or not the CME current undergoes 
a non-perturbative renormalization, or whether it survives in the strong cou- 
pling limit, is still an open question [^Olf¥T] . We will characterize the real-time 
dynamics of topological charge in the system by a certain distribution in the 
chiral chemical potential ^5. 

This paper is organized as follows. In Sec. [2] we briefly review the experimen- 
tal observable accessible at present in heavy ion collisions and motivate the 
connection between this observable and the electric-current correlation func- 
tion and susceptibility. The following Sec. [3] is divided into four subsections: 
first, we re-derive our previous result for the CME current [3T] and compute 
the longitudinal susceptibility using the thermodynamic potential. Calculating 
the transverse susceptibility is not so straightforward as the longitudinal one, 
and we explain our method of computation in the three following subsections. 
We then proceed to Sec. H] where we make a comparison of our results to the 
lattice QCD data and provide the formulae that can be used in the description 
of experimental data. Finally, we summarize in Sec. 



2 Experimental Observables and Correlation Functions 

In heavy ion collisions at any finite impact parameter the distribution of mea- 
sured (charged) particles depends on the azimuthal angle 0. The hadron mul- 
tiplicity as a function of can be decomposed into Fourier harmonics: 




oc 1 + 2v 1± cos(A0) + 2a± sin(A0) + 2v 2 ± cos(2A0) H , (1) 



where N + and iV_ are the number of positively and negatively charged parti- 
cles respectively, and A0 = — \1/ R p is the angle relative to the reaction plane. 
The coefficients V\ and v 2 quantify the strength of "directed" and "elliptic" 
flows respectively; they are not expected to depend on whether measured 
particles are positively or negatively charged. On the other hand, the term 
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proportional to a± is V- and CV- odd and describes the charge separation 
relative to the reaction plane - in other words, as sketched in Fig. [TJ a± quan- 
tifies the strength of the charge flow directed perpendicular to the reaction 
plane. Unlike v\ and t>2, as we will see later, a + and a_ depend on the charge 
carried by measured particles. If the CME electric current is directed upward 
(downward), we expect a + > (a_ < 0); non-zero values of a + and a_ indicate 
the presence of V- and CP-odd effects. 





\ Vi / 









Fig. 1. Collision geometry and collective flows decomposed in Eq. ([T]), of which a± 
is sensitive to V- and CP-odd effects. 

We can carry out the decomposition ([1]) for each event and take the ensem- 
ble average over all the events. We shall denote this averaging procedure by 
((■ • ■)) throughout this paper. Because the topological excitations fluctuate not 
only locally (point-by-point in space) but also globally (event-by-event), ((a±)) 
becomes zero and the symmetry is restored in a sense of average. It is nec- 
essary, therefore, to measure the correlation functions ((a + a + )), ((a + a_)), and 
((a_a_)), which are invariant under V and CV transformations and thus their 
ensemble average is non-vanishing. 

The experimental observable sensitive to ((a±a±)) was proposed by Voloshin [32] 

1 



«cos(A0 a + A(j)p))) = 



cos(A0 c 




Here a and (3 indicate either + or — charge of measured particles, and the 
sum goes over all charged hadrons in a given event. This observable has an 
important property that becomes clear when one explicitly isolates the terms 
B a p driven by fluctuating backgrounds 



(3) 



((cos(A0 a + A(j)p))) = ((cos A(f> a cos A<fip)) — ((sin A0 a sin A(pp)) 
= (((vi^)) + B%) - (((a a ap)) + B$) . 



If the in-plane B 1 ^ and out-of-plane 5°^ backgrounds are the same, they 
cancel out and since for a symmetric heavy ion collision in a symmetric rapidity 
cut ((vi !0c vi t p)) ~ 0, one can identify the above (j3J) with —{{aadp)) that serves 
as an order parameter for V- and CP-odd effects. 

Let us now take a further step to establish a link between the above experi- 
mental observable and the theoretical computation. In theory a clear manifes- 
tation of the CME is spontaneous generation of an electric current, J, along 
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the direction parallel to the applied magnetic field, B [30f3l] . We can antici- 
pate that this electric current of quarks, if J z > 0, would induce a + > and 
a_ < in Eq. ([T]) for the observed hadron distributions. This is due to i) the 
quark-hadron duality and ii) the rapid transverse expansion of the system that 
prevents the charge asymmetry from being smeared out by thermal diffusion. 
It seems natural that a theoretical quantity relevant to ((a±a±)) should be a 
correlation function or fluctuation with respect to the electric current. For one 
event it would be legitimate to accept the following relations: 





N- 




^cosA0 +ii « 




OC J x 


t=i 


1=1 




N+ 


N- 




sin A0 + j « 


Sil1 A< ^-,i 


OC J 2 


3=1 


i=l 





Indeed, J2 cos A0 +i ;, f° r example, counts an excess of positively charged parti- 
cles in the region near A<ft « (i.e. moving in the x > direction) as compared 
to that near A0 w 7r (i.e. moving in the x < direction). Here we defined a; 
and z as the transverse and longitudinal coordinates, respectively. From now 
on, to make our notation more general, we will use J± and J\\ instead of spe- 
cific coordinate indices; namely J± is a current component perpendicular to 
B and J\\ parallel to B. We can then identify 

cos(A0 a + A^)oc-^-(ji-Jf). (6) 

Note that we have neglected the screening (absorption) effect by the medium 
for the moment; we will discuss the medium screening in more detail in Sec. HI 

To connect these expressions to the experimental observables we should further 
take the average ((• • •)). Computing this event average would however require 
a complete real-time simulation describing the dynamics of topological transi- 
tions in the produced matter. Since such a simulation goes beyond the scope 
of the present work, here we will adopt the following working definition for 
the correlation function that is related to the experimental data. 

(1) First, we identify the fluctuation measured on the event-by-event basis 
with the spatial correlation function obtained in the field theory calcula- 
tion. This is a common assumption used for example also in discussions 
on the QCD critical point search (56]. If the system's volume is sufficiently 
larger than the correlation length, this assumption should make sense. 

(2) Second, we compute the correlation function for a fixed value of the chi- 
ral chemical potential ^5 that induces a difference in number between 
right-handed and left-handed particles; N 5 = N R — N L . 
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(3) Third, we make a convolution of the results at fixed /is with a weight 
function of the /is distribution. We will adopt a Gaussian ansatz for the 
weight function; the dispersion is characterized by the rate of topological 
transitions depending on the system's temperature. 

One may wonder whether the above-mentioned procedure is a correct one, or 
merely a convenient assumption. Suppose that (J]_) and (J|) were computed 
from the first principles of quantum field theory avoiding the steps (2) and 
(3) above, as is the case for the lattice QCD simulations. Then one might 
think that the results are to be interpreted as quantities directly relevant to 
the experimental data. We think however that the situation is not so sim- 
ple. The difficulty lies in the different dynamics of topological excitations in 
Euclidean and Minkowski space-times. In the case of Euclidean space-time 
the instantons are exponentially suppressed by the Boltzmann factor, while 
in Minkowski space-time the sphalerons are not [57]. Therefore, it is essential 
to take account of sphaleron-like effects in real-time dynamics. If we make 
use of the imaginary-time formalism of the finite-temperature field theory, it 
would describe properly the instanton transitions but would miss the unsup- 
pressed sphaleron contributions unless the analytical continuation is correctly 
formulated. Usually the analytical continuation from the imaginary- to the 
real-time dynamics is not obvious though. It would thus be appropriate to 
solve the real-time dynamics e.g. by using the spatially discretized Hamilto- 
nian [58,59, 6T)|61ll(j2j in the presence of external magnetic field. This is how- 
ever quite difficult technically and hence, in this work, we have employed a 
pragmatic approach assuming that the computation of the spatial correlation 
function at finite /A5 and the computation of the average over the distribu- 
tion in fi 5 (that is governed by the real-time dynamics) are separable. Such 
a separability would hold for example if the CME response of the system to 
a finite /is were much faster than the dynamics responsible for the genera- 
tion of the distribution in /is, and the back- reaction of the CME current on 
topological transitions could be neglected. We think that this is a plausible 
assumption since the generation of the CME current in response to a finite /is 
is not damped by interactions with gluons, and at least at large N c the back- 
reaction of the CME current on the topological dynamics of gluon fields can 
be neglected. In this approach, the properties of the topological transitions 
in hot QCD medium such as the sphaleron rate are incorporated through the 
weight function for the /is-distribution. 

The correlation function, in general, has two contributions: one is from a 
disconnected (tadpole-tadpole type) diagram and the other from a connected 
diagram. The former is related to the expectation value of the induced current, 
and the latter is given by the corresponding susceptibility. Here we introduce 
the following notations: 

(JiU = (J\\)l 5 + xl , (Ji), 5 = (J±)l 5 + xi ■ (7) 
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We note that (• • -) M5 and above are an expectation value and suscepti- 
bility of the electric current, respectively, for a fixed value of /15. Under our 
assumptions listed above, the ensemble average reads, 

((cos(A0 + + A0 + )» = -((cos(A0 + + A0_)» 
= ^/d/i 5 %)((^-(^) 

= ~M I d/i5 w M< J n>« + xl - xi) , (8) 

where W(/xs) is the weight function normalized such that / dfi 5 W(/x 5 ) = 1 
(and thus W(/is) has a mass dimension in our definition). We will specify this 
function later. The overall constant, c, is treated as a free parameter which is 
independent of collision geometry, and we approximate N + »s AL « JV±. 

It is interesting to see that the combination of the transverse and longitudinal 
current correlators in Eq. (jSj) has a clear correspondence to the form of Eq. (j3J). 
In (Jj)^ — (Jf)n 5 together with Eq. (jjTj), there is a cancellation between major 
components in and xji 5 > that is interpreted as the cancellation between 
backgrounds B™a and B^g discussed above. However, because of the presence 
of the external magnetic field the cancellation is not necessarily exact, which 
gives a background on top of the CME contribution from (</||)L- 



3 Current and Susceptibility 



If we neglect the quantum fluctuation of gauge fields on top of the external 
magnetic field B, we can explicitly accomplish the integration in the quark 
sector with full inclusion of B. We can write down the thermodynamic grand 
potential at finite T and fi q , which is defined by Q = —T In Z from the partition 
function Z, given as 



/ 



a: 



ks 



s=± k=0 



00 dih r 



^fcA+^5Z ln ( 1+e 



, (9) 



where /, s, and k refer to the flavor, spin, and Landau level indices, respec- 
tively. The factor a{ s is defined by 



for k = and qjB > 
for k = and qjB < 
for k > , 



(10) 
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which takes care of the spin degeneracy depending on the zero or non-zero 
modes. The quasi-particle dispersion relation is 



u kX = ^(y/pi + 2\qfB\k + \ f x^ + M|. (11) 

Here A = sgn(p 2 )s stands for the helicity and My represents the quark mass for 
flavor /. We note that u>k\ is a flavor dependent quantity, though we omit an 
index / for concise notation. In our convention p z refers to the z-component of 
the three momentum p (not the third component of p^ which has an additional 
minus sign from the metric). 



3. 1 Differentiating the grand potential 



We can calculate the expectation value and the susceptibility of the electric 
current operator by taking a functional derivative with respect to the (source) 
gauge field as discussed in Ref. [31J. In general the electric current can be 
expressed as 

ST[A] 



0" M (*)> 



(12) 



A=A 



where T[A] is the effective action in a certain gauge and A represents the 
background gauge field (corresponding to the external magnetic field in our 
case). In the finite-temperature field theory the effective action (potential) 
translates into the thermodynamic potential. Thus, the induced current is, 



<J") 



, 5T\nZ[A] 



d 4 x 



SA^x) 



A=A 



d 4 x 



SQ[A] 



SA^(x) 



(13) 



A=A 



It is easy to write down an expression for the current in configuration space 
using the above equation with the functional derivative, which turns out to be 
equivalent with the diagrammatic method (see Appendix |A|). Because (j M (x)) 
does not have x dependence under a spatially and temporally homogeneous 
magnetic field, we can replace the functional derivative by a derivative with 
respect to homogeneous A^, that is, 



<J") 



dQ[A] 



dA, 



(14) 



A=A 



The advantage of this rewriting is that we can now directly work in momentum 
space. Then, because of the structure of the covariant derivative p—qfA (where 
qf > for positively charged flavor), the derivative with respect to A^ can be 
replaced by that with respect by p^ times —qf. To evaluate the current we 
have to take the derivative with respect to A 3 which becomes the derivative 
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with respect to — p z times —qj. In this way the current parallel to the external 
magnetic field B is immediately written down as [31] 



Here we have used the fact that the spin sum (s = ±) makes a cancellation for 
non-zero modes and only the zero-mode contribution from k = remains non- 
vanishing. The current has an origin in the quantum anomaly coming from 
the surface term of the p z -integration, so it is an exact result and insensitive 
to any infrared scales such as the temperature T, chemical potential fi q , and 
quark masses Mf (see also Refs. [63|64II65] ). Since rotational symmetry in the 
transverse plane perpendicular to the magnetic field is kept unbroken, the 
transverse currents are zero; (J±) fl5 = 0. What we have addressed so far is 
just to remind the discussions given in our previous work [31] . We present an 
alternative derivation of the induced current using the exact propagator in a 
magnetic field in Appendix [A] 

In contrast to the current expectation value, the (unrenormalized) current 
susceptibility is dominated mostly by non-zero modes at high Landau levels; 
this makes the susceptibility (but not the difference of parallel and transverse 
susceptibilities) sensitive to the UV regularization, making it unphysical. In- 
terestingly enough, the zero-mode again plays an essential role in a finite 
difference of x ~ X~^ 5 ■ The susceptibility is deduced from 

* = -T™ (16) 

X g A 2 • 

in the same way as the previous discussions on the induced current. The 
parallel component is again easy to evaluate because the structure of the 
covariant derivative is unchanged. We must here introduce a UV cutoff A. 
We shall impose A in a symmetric way, that is, by Heaviside's step function 
6>(A 2 — p 2 z — 2\qfB\k). We note that our final result of x ~ xjl 5 is UV finite and 
thus it does not depend on the cutoff parameter nor scheme. Our naive way 
to cut the momentum integration off is not gauge invariant, and so it is of no 
use in order to compute x and x^, individually. 



In our prescription the ^-integration is bounded by A& = yA 2 — 2\qfB\k 
for a given k, where k takes a value from to k\ = [A 2 / (2\qfB\)\ . Thus, 
by differentiating the grand potential with respect to A z twice, we have the 
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longitudinal current susceptibility as 

/ f Ak d P* d2 r lT n (-\ j_ - 



X 



K a ±m,)/T 



,21 



^ Qf\qfB\ f p z A/i 5 \r _ 



/-<;- 



V7W C £ + £M [, _ „ fK ,) _ S>A>) ] , (17) 

/,»,* ^AA A / L 



where cja s = y (A + S/U 5 ) 2 + M| and the Fermi-Dirac distribution functions 

arera F (co') = [e^-toW 1 " and n f (u;) = [e^ u+flq ^ T . We have introduced 
a new notation which is the spin degeneracy defined by 



9k = + oiks 



1 for k = 

2 for k ^ 



In principle, the current susceptibility has an additional contribution from 
mixing with the chiral susceptibility xm through dependence of the dynami- 
cal mass Mf on A^, resulting in a contribution like (dMf /dA^) 2 XM- This is, 
however, negligible for the strength of magnetic fields and associated currents 
relevant to heavy-ion collisions. (We have numerically confirmed this by using 
a mean-field chiral model.) 

The longitudinal susceptibility (|17p is UV divergent, as we have already men- 
tioned, and is strongly dependent on the value of cut-off A and how A is 
imposed. In fact it is straightforward to take the limit of B = n§ = in 
Eq. (|T7|) . which leads to an orientation independent susceptibility; 

,2 A3 



VTN C J2 I 1 ~ n r ^ " n F M] > ( 19 ) 



Xo ■ - -I- / , 

j 67T Z UJ\ 



where we have defined u\ = J A 2 + M 2 . This is UV divergent proportional 
to A 2 for large A. It is necessary, therefore, to perform the renormalization 
to extract a finite answer for x — Xo-, i n which the leading divergence ~ A 
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cancels but the logarithmic divergence may remain in general. As long as xL ~ 
X^ 5 is concerned, as we are explicitly computing, only a finite term appears 
and there is no subtlety associated with UV divergence and renormalization. 

Evaluating Xu 5 is n °t as easy as xi s ] we cannot simply replace p by one aug- 
mented with A in the transverse direction because of the Landau quantization. 
In what follows below we will explain how to calculate this quantity by taking 
three steps. 
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Step i) We should return to the original definition of the current correlation 
function in terms of quark fields; Xi °t Jd^xd^y^Y^i^^Y^iv)) ■ We already 
know the answer when B = fi 5 = 0, which is given in Eq. f)19p . In the first 
step, hence, we will confirm this known answer by the diagrammatic method. 

Step ii) We will introduce an external magnetic field B in the second step. 
The quark propagator requires a modification with the projection operator in 
Dirac indices. Then we will find that the transverse susceptibility takes the 
same form as the longitudinal one except for the Landau zero-mode. 

Step iii) Finally we will extend the calculation to the case with finite To 
compute the quark propagator we need to insert one more projection operator 
with respect to the Dirac indices, which separates opposite helicity states. 

3.2 Establishing the diagrammatic method — step i) 

In Minkowskian convention we can express the electric-current correlation 
function as 



By taking the contraction of quark fields, the above expression is decomposed 
into the disconnected (tadpole-type) and connected contributions. 



Fig. 2. Diagram of the connected contribution to the electric-current correlation 
function (photon self-energy) and the momenta running in the loop. 

The diagrammatic representation of the connected contribution from Eq. (f20l) 
is displayed in Fig. [2j from which we can recognize that this is nothing but 
the diagram for the photon self-energy (polarization tensor, see Ref. [55] for 
discussions on the vacuum polarization in a magnetic field). Therefore the 
answer must be zero in the limit of zero momentum insertion from the external 
legs; otherwise the photon becomes massive and breaks gauge invariance. It is 
obvious from such a recognition that the UV divergent answer in Eq. ffl9|) is 
unphysical and merely an artifact from the naive cutoff prescription. 

The corresponding expression to Fig. [2] translates into the current susceptibil- 




(20) 




P = (Po,p) 
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ity of our interest here given in the following form; 



(21) 



where G(x,y) represents the quark propagator in configuration space. Inte- 
grating over x and y and using the familiar expression of the propagator in 
momentum space for B = /15 = (then there is no preferred direction and so 
we put a subscript instead of i) we immediately reach 



x ^ 2 f d 3 p f T d Po 



(2tt)' 



tr 



7 



— Mf p 1 — Mf 



v 2 f d 3 p [Tdpo 4(2p iqi +p-q-M 2 f ) 

1V c Y f J 2vr (pS-^XpS-^) 



(22) 



where p = (p = 2ni(n + |)T + /i g ,p) and q = (p , q — > p) and / (dp /2 
iTJ2 n - We can carry out the Matsubara frequency sum over n to find, 

. f Tdp 4(2p iqi +p-q-M 2 f ) 
2vr (p§-u;;)(p§-a£) 

1 ^g(^ F (^)+^ F (^))-^p(^ F (^p)+^ F H)) 

2(u; p + cj 9 ) u; g - u; p 

(2^% -p-q-Mj) [u- 1 (n F (u q ) +n F (w g )) - w~ 1 (n F + 



q^p 



2(U p + UJg)(U g - Up) 

i 2 d 



,_d_ \Pi_ 

dpi \ io v 



1 - n F (uj p ) - n F (ujp) 
1 - n F (u)p) - n F (u)p) 



+ 



Pi 



1 - n F (ujp) - n F (u p ) 



ut 



du r 



(23) 



Plugging the above (1231) into Eq. (122]) and integrating over the momentum p 
with a UV cutoff A, we can correctly rederive the result of Eq. ( fl9l) in the 
diagrammatic way. 



Before closing this subsection we should mention on our procedure to take 
the limit of q —>■ p. In this subsection this procedure is not indispensable; we 
could have started the calculation with q = p. Then, the residue of a double 
pole gives a term involving the derivative of the distribution functions with 
respect to u p which is exactly the same as the term arising from the q — > p 
limit. The reason why we distinguished q from p is that, as we will see shortly, 
the presence of B ^ induces a difference in the transverse component in 
Landau-quantized momenta. Thus, the calculation process elucidated in this 
subsection has a smooth connection to the later generalization. 
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3.3 Turning on B — step ii) 



In the presence of non-zero B the transverse and longitudinal directions be- 
come distinct even before introducing finite ^5, which should result in Xo ~ 
Xq 7^ regardless of topological excitations. We already know the answer for 
the longitudinal susceptibility from Eq. (fTTI) . that is, by taking the limit of 
/is = 0, we can readily write down, 



}% 2n 



1 -n F (uo k ) -n F (u K ) . (24) 



Let us calculate the transverse susceptibility by means of the diagrammatic 
method developed in the previous subsection. For this purpose we should first 
come by the quark propagator under a constant magnetic field. 

The solution of the Dirac equation with a constant magnetic field is known, 
which forms the complete set of orthogonal wave-functions. We here introduce 
several new notations. We use the following basis functions in a gauge choice 
where A = A x = A z = and A y = Bx; 

f k+ (x) = <j> h {x- Py /(qB)) (k = 0, 1, 2, . . . ) , 

f k .{x) = <f> k . 1 {x-p v /{qB)) (A; = 1,2,3,...). 1 ] 

with 4>k(x) being the standard Landau-quantized wave-function defined by 

= \f^n(—) ' exp(-l-\ q B\x 2 )H k (J\a~B\x) , (26) 



, 2 k k\ V vr J "I 2 

where H k (x) represents the Hermite polynomial of degree k. (We omit the sub- 
script / of qf for the moment.) We can easily confirm that the basis functions 
satisfy the orthogonality property as follows; 

dx f k+ (x)f l+ (x) = 8 K i , 
dxf k+ (x)fi_(x) = (I > 1) , 

(27) 

dxf k _(x)fi_(x) = 5 kj i , 
dx f k _(x)fi+(x) = 4-1,2 {k > 1) ■ 

We can define the projection matrix with respect to the Dirac index according 
to Ritus' method [671 as 



f k+ (x) + f k 4x)\+-[f k+ (x)-f k 4x)\ 1 1 7 2 J (28) 



Pk{x) = \ 
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for qB > and f k+ and f k - are swapped for the case with qB < 0. It is then 
straightforward to prove that 



(id-qA- Mf)P k (x) e -' l{pot - p y y - pzZ) 

(29) 

= P k (x)(p ol ° + sgn^^lqBlk^ -p zl 3 - M f ) e -^t-Pyy-p^) ; 
where the right-hand side is expressed in a form of the free Dirac operator with 



a modified momentum p = (po,0 } —sgn(qB)y2\qB\k,p z ). Hence, the solution 
of the Dirac equation with a magnetic field is given by a combination of the 
projection matrix and the free Dirac spinors u(p,s) for particles and v(p, s) 
for antiparticles with a momentum argument p, that is, 

P k (x) u{p, s) e -KP0t-p y y-Pzz) ^ v ^ g j e i(pot- Py y- Pz z) _ ( 30 ) 

Here we note that P k (x) is a real function. Using these complete-set functions 
and defining associated creation and annihilation operators, we can compute 
the propagator (ip(x)ip(y)) to reach 

G(x,y) = {ip(x)i[)(y)) 

dpo r dp y dp z e -i\po(xo-yo)-py(xy-y y )-p z (x z -y z )] , 31 x 

2tt J (27r) 2 - 1 

xP^i^-M/)- 1 ^). 

It should be noted that P k (x) and P k (y) have implicit p y dependence in the 
argument (see Eq. (|25|) ). The spin degeneracy is automatically taken into 
account by the property that fo-(x) = and thus only one spin state has a 
non-zero contribution for k = 0. 

We use this form of the quark propagator to express the susceptibility (|2~T|) . 
Then we have to perform the coordinate integration with respect to xo, yo, 
x x , y x , x y , y y , x z , y z , and the momentum integration over p , p' , p y , p' y , p z , 
p' z and also the summation over k and /, where p' , p' y , p' z , and / are from 
another G(y,x). As usual, after we integrate with respect to x , x y , and x z , 
we have (2n) 3 S(p — p' )S(p y — p' y )6(p z — p' z ) (there is no momentum insertion 
from external legs), from which only the integrations with respect to po, p y , 
and p z remain. Then, renaming x x — > x and y x — > y, we have 



x tr 



fP k (x)(p - MfY'PM-fPtiy)^ - MfY'P^x 



where we have defined q = (po, 0, — sgn(qB)y2\qB\l ,p z ), L y = J dy y , and 
L z = j dy z . Because p y is now common in all P k (x), Pi(y), etc, a shift in the 
integration variables, x and y, can get rid of p y from the integrand. Therefore, 
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J(dp y /27i) simply counts the number of Landau-quantized states, that is given 
by L x (\qB\/2ir). Here L x L y L x is nothing but the volume V. After all we can 
write the above into a form of 



Xi 



iVTN c J2 



q 2 MfB\ 



2tt 



2tt J 2tt 



dx dy 



x tr 



yP k (x)(p- M f Y 1 P k (y)^ Pi(y)(^ - M f Y l Pi 



(33) 



Apart from the projection operator which we will discuss soon below, this 
expression has a diagrammatic representation with momenta running in the 
loop as depicted in Fig. [31 



P = (Po, 0, -8ga(qB)y/2\qB\k,p z ) 




Fig. 3. Diagram and the running momenta corresponding to Eq. (|33p . The Landau 
levels denoted by k and I are constrained (either A: = /or/c = / + lor/ = A; + l)by 
the projection operators. 



Because Pk{x) is a combination of the unit matrix and 7 X 7 2 as defined in 
Eq. (1281) and thus is commutative with 7 3 , the longitudinal susceptibility is 
easy to evaluate. Using Pi(x)j 3 Pk(x) = Pi(x)Pk(x)-y 3 and the orthogonal- 
ity ( 1271) . after integration over x, we obtain 5i t k for k > and (1+i sgn(g5)7 1 7 2 ) /2 
for k — I — 0. The calculation is then reduced to the same as in the previ- 
ous case with several replacements; the momentum p by p, and the integration 
with respect to transverse momenta by the summation over the Landau levels. 
Then the longitudinal result of Eq. ( 12^1) is almost trivially understood from 
Eq. (T23D with i = z. 



We are now ready to consider the transverse polarization. In this case 7 1 and 
7 2 are not commutative with Pk{x). Because the calculations for 7 1 and 7 2 
are just parallel, we focus only on the case of 7 1 below. Using 



dxP^x^Pkix) 



( 1 



x 1 1 x 2 

-0l-i,kl - ~0l-i,kl 

1 r 1 1 r 2 

-<)z,fc_i7 + -fy,fc-i7 



I 2 



$l-l,k + Olk-l 



7 



8l-l,k — 5; fc_i 



for k = 
for / = 
7 2 for fc 7^ 0, / ^ 



(34) 
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we can evaluate the trace with respect to the Dirac index to arrive finally at 



x ^ = WTN c y^^Y, — — , \77 7 ■ (35) 



Up j 

It should be mentioned that there emerges no spin degeneracy factor in the 
above expression. Taking the Matsubara sum leads to 

Tdpo 2{p 2 -p 2 z -M 2 )5 l ^ k 

2vrV (p8-w?)(p8-^) 
A; + 1 r i k 



UJ. 



1 & 

- 1 - 72 F (^) - n F (^g) 1 - 71,. (Wp) - ^ F (^p) (36) 



g - -Up 



with, as indicated in Fig. [31 p = (po = 2m(n-\-\)T-\-[i qi 0, — sgn(qB)J2\qB\k,p 



and g = (p , 0, — sgn(qB)^j2\qB\l ,p z ), where I — k + 1. Here we note that 
and I are both bounded by the UV cutoff, that is, the maximum k is k — 1 
with = [(A 2 — pl)/(2\qfB\)\, and p 2 is cut off by A. We realize then that 
subsequent terms are canceled in summation over k and only the edge term 
(from k = k — 1) remains non- vanishing, which is reminiscent of the remaining 
surface term in the B = case out of the momentum integration of Eq. ( j23|) . 
Therefore, at last, we get 



Vo = VTN C Y: r dp, — [l - n F (u A ) - n F (co A )} (37) 

- Z7T J— A id A 



after the summation with respect to k from to k — 1. To be more precise, 
it should be mentioned here that we approximated as 2\qjB\k + p 2 z + M| — > 
u\ = A 2 + M 2 , which is valid as long as A 2 3> \qjB\. 

It is very interesting that we can find a simple analytical formula for the 
difference Xo ~ Xo- From Eqs. fl2M and fl37|) we can write, 



Xo ~ Xo 

vr^E^-fi-",K)-fl>A)i^-«- -' ' A - ^ (s 




where, as we have already defined, A k = y A 2 — 2\q$B\k and k = [(A 2 — 
p 2 ,)/(2\qfB\)\ . The point is that we can rewrite the latter term in the paren- 
theses into a form of the discrete summation from the fact that k is defined 
by the floor function. That is, 

A _ rA rA 1 rA 2 

dp z k = dp z x + / dp z x 1 + / dp z x 2 H 

JAi JA 2 JA 3 /ggN 

(A, - A 2 ) + 2(A 2 - A 3 ) + 3(A 3 - A 4 ) + • • • = £ A k . 

k=l 
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This form (times two) exactly looks like the first term in the parentheses except 
for the k = term. Therefore, only the zero-mode with k = contributes to 
the final expression; 

Xo ~ Xo = VTN C qf ^ f ^ — i 1 - n F M - n F (w A ) 

(40) 

Although Xo an d Xo are divergent ~ A 2 in the gauge-variant cutoff scheme, as 
we have mentioned before, the difference between them is certainly UV finite, 
so that it can be a well-defined quantity. 



3.4 Introducing /z 5 — step Hi) 



To treat the situation in the presence of not only B but also ^5, we need to 
insert the projection matrices defined as 



r ± (p) = i(l±p- 77 7 5 ) 



(41) 



to calculate the quark propagator which involves the inversion of Dirac ma- 
trices. The above-defined projection matrices have the following property, 



- U> + /"57V + M)(p 2 - 2 M5 p-77°7 5 - £ - M 2 )~ l T ± (p) 



(^ + /i57°7 5 + M) pl-{\p\±^) 



r±(p). 



(42) 



We can identify ± as A appearing in Eq. (j5J), so that T±(p) is the projection 
operator to the state with helicity A = ±. We insert the unity, 1 = T + (p) + 
r_(p), after the propagator to take the inversion. By doing this we can express 
the susceptibility in the following way; 



Xi 



iVTN c J2 



2tt 



E 

k,l 



dp z r T dp 



2tt 



2tt 



dx dy 



x tr 



fP k (x)(p + /i 5 7V - M f )- L (T + (p) + T4p))P k (y) 
x fPi(y)d + /i 5 7°7 5 - M f )-\T + (q) + r_(g))P,(x) 



(43) 



Since we know the answer for the longitudinal susceptibility, let us concentrate 
in calculating the transverse one (i.e. i — x) hereafter. After the integration 
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over x and y with Eq. ([54"]) we can find 



A 



iVTN c J2 



2tt 



E 

A-./ 



X 



E 

A,A'=± 



2tt 7 2tt 
^a A , + T||, + iT A |, - iTf A , 



(44) 



where 



1 w — u 



YiP + /^57°7 5 + M/)r A (p)y (| + /i 5 7 U 7 5 + M f )r A ,(g) 



0„,5 



(45) 



We note that tr[7 M 7 v 7 p 7°"7 5 ] = — 4ie Ml/pcr in the convention we are using (where 
+1). After some lengthy calculations we can notice that T AA , becomes 



as simple as 

= {d ij + AA'[(p l g j + fPq 1 ) - 5 i] p ■ q] + ie 0i3k {Xp k - A'g fc )} 



x 



p 2 + XX'(\p\ + Xfi 5 )(\q\ + X'fis) - M 



(46) 



from which we get 



rpll | rp22 , -T-il2 -rp21 



pi + AA'(|p| + X^)(\q\ + AV 5 ) - (1 - Ap 3 )(l + A'g 



The Matsubara sum amounts to 

" ' 2, 2 (p8-wjj(pg- w ; A o 



(47) 



A,A'=± 

E 



a,a'=± 2u; gA '(^A' ~ u px 



|5| + A'/is) 2 + AA'(|p| + Xfjt s )(\q\ + X%) 



x (1 - Xp s )(l + A'g 3 ) 1 - n F (u q y) - n F (oj qX >) 



+ E 



X,X'=± 2w pA(^pA - ^gA'/ 



\p\ + A/i 5 ) 2 + AA'(|p| + Xfi 5 )(\q\ + X'fj, 5 ) 



x (1 - Xp s )(l + A'g 3 ) 1 - n F (oo pX ) - n F (w pA ) 



(4* 



Here we note that this expression consists of two parts; the former part 
contains the distribution function with an argument u q x> (i.e. n F {uj q x>) and 
n F {uj q x')) and so the summation over A can be easily taken to simplify the 
coefficient in front of [1— n F {oj q x>)— n F (uj q x>)]- Remarkable simplification occurs 
as a result of the helicity sum, which leads us to 



Eq. (BED = £ ^ 



1 + 



1 - n F (uj q x>) - n F (uj qX >) 



E — ( 1 + 7=!?) f 1 ~ M^pa) - M w p a) 



A=± 



<^pA V 



\P\ 



(49) 
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From this it is apparent that there is a major cancellation in the summation 
with respect to k and only the edge terms remain non-vanishing. Hence, we 
carry the summation out to find the following expression; 

^E^TT / A dp z — (l + S ^)\l-n F (uj AX )-n F (uJAx)} . (50) 
77 An z J -a ujax V A / l j 



Therefore, after lengthy procedures in this subsection, what we find out at 
last is almost the same as Eq. (140j) with a minor modification by /is, that is, 



2U o, (51) 



A-oo y 27T 2 

In the limit of A — > oo the susceptibility difference has no dependence on ^5. 



5.5 Discussion 



Since the final results are so simple, we will present another way to get the 
same answer from heuristic arguments. Let us consider a situation in which a 
homogeneous magnetic field B is parallel to a homogeneous electric field E. 
In such situation a current parallel to the magnetic field will be generated by 
the Schwinger process. For infinitesimal E = (0, 0, E z = E) the rate of change 
of this current is determined by the electromagnetic anomaly relation leading 
to 

(g)-™.£*i!g!* +0 c*>. 

We can choose a gauge so that E = — 8qA z , and then A z is also infinitesimally 
small, which allows us to use the linear response relation to express the current 
changing rate in another form, 



- / d 3 xd V/^^ i||(a /)\ A z (x') + 0(A 2 z ). (53) 



\dx / J \ dx I re t 

Assuming the translational invariance in time in (jn(x)ju(x')), we can replace 
d/dxo acting on it by — d/dx' , and then we can perform the integration by 
parts to move d/dx' acting onto A z (x') which results in E = —8qA z . Eventu- 
ally we arrive at 

A T \ 

d 3 x dV OH (x) n (x% et E + 0(E 2 ) . (54) 
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Identification of the left-hand side of Eq. (1521) with that of Eq. f|5^|) concludes, 



= VTNc ^_ME\ (55) 



Now we could do the same derivation for ( J]_) to find that it is zero. The rea- 
son is that (dj±/dxo) is always vanishing for infinitesimal E in the transverse 
direction; small transverse electric field cannot give rise to an electric current 
in the transverse direction because of the energy barrier by the Landau quan- 
tization. Therefore x -1 — 0, and so Eq. (1551) gives the difference x ~ X~S which 
is in agreement with our results (|40l) or (}5Tj) . 

The anomaly equation is an exact relation. Therefore, since this derivation 
clearly shows that the electric-current susceptibility is determined by the elec- 
tromagnetic anomaly, we conjecture that at least for massless quarks our re- 
sults are exact and will not be modified by including perturbative gluonic 
interactions that do not affect chirality. 

In Appendix [B] we will argue that the current susceptibility is in some sense 
similar to the current-chirality correlation. There, we will see that the anomaly 
relation again constrains the current-chirality correlation. 

Our results and the above-mentioned heuristic arguments are consistent with 
what is known in condensed matter physics. The hall conductivity (the 
electric conductivity in the z-direction when E is imposed in the j-direction) 
in the presence of the magnetic field is a quite familiar quantity and our x 
and X" 1 could be regarded as a zz and a xx in the terminology of the quantum 
hall effect apart from the fact that the susceptibility is a quantity in the zero- 
momentum limit at zero frequency, while the conductivity can be a function 
of finite momentum and frequency [54J . Because our calculation has three spa- 
tial dimensions, the counterpart in condensed matter physics is the multilayer 
Dirac electron system. Then, as discussed and confirmed in Ref. [68], only the 
Landau zero-mode dominates the physics properties; a xx is zero because the 
transverse current operator involves a shift in the Landau level as is embodied 
in Eq. in our calculation, which is actually a common knowledge in con- 
densed matter physics. The longitudinal one, o zz on the other hand, is finite 
and is basically given by transport from zero-mode on one layer to zero-mode 
on another layer. In our calculation for relativistic quark matter it is intrigu- 
ing that such a finite contribution along the longitudinal direction is uniquely 
constrained through the exact anomaly relation. 
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4 Discussions on Lattice QCD and Experimental Data 



Now that we have a simple expression for Xo — Xo > h would be an interest- 
ing question whether our results are consistent with the existing results from 
the lattice QCD simulation [37] and from the previous work [30] aimed at 
describing the CME in heavy ion collisions. In this section we will see that 
our estimate shows reasonable agreement with lattice QCD and the formulas 
given in the previous work. 



4-1 Lattice QCD Data 



In Ref. |37J the ITEP lattice group reported on the electric-current susceptibil- 
ity in quenched QCD which they found to be larger in the longitudinal direc- 
tion than in the transverse direction. The simulation condition is that Nf = 1 
with q = — e/3, the lattice spacing a = 0.095 fm (i.e. A = n/a = 6.51 GeV), 
the volume L = 16a = 1.52 fm (and thus V = 460 GeV -3 ), and the simulation 
is done with N c = 2. Our calculations make sense above T c where the physical 
degrees of freedom are quarks, so let us compare with the lattice results at 
T = 1.12T C = 350 MeV. 



First of all, plugging the above numbers into Eq. (1T91 we find Xo as 
Xo 2TA 2 2 x 0.35 x 6.51 2 



q 2 V 2 3n 2 V " 3tt 2 x 460 



GeV b = 2.2 x GeV b , (56) 



which is much larger than the lattice result. This is not contradictory, however. 
As we have emphasized repeatedly in this paper, a non-zero value of Xo is an 
artifact of the naive momentum cutoff, and our scheme of UV regularization 
is different from the lattice one. 

More interesting is the comparison of the U V-finite difference xl ~ Xo > which 
turns out to be 

= = 2 X °f * GeV 4 = 7.7 x 10-*|,B| GeV 4 . (57) 

q 2 V 2 2n 2 V 2tt 2 x 460 |y 1 V ' 

This is fairly close to the lattice results of Ref. [37] ■ To make the comparison 
easily visible we shall compare our estimate with the lattice data corresponding 
to Fig. 8 in Ref. |37J for the difference between the longitudinal and transverse 
susceptibilities. Figure H] plots our estimate of Eq. (1571) versus the lattice QCD 
data. It is apparent that Eq. (1571) is in good agreement within the error bars 
of the lattice data. 

One may wonder why Eq. (J57J) works so well even though it does not contain 
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Fig. 4. Comparison between the lattice QCD data from Ref. [37] (shown by the 
square dots) and our estimate (shown by the solid line). 

any information on the topological excitations and thus no dependence on ^5. 
We think that a possible explanation is the following: Euclidean lattice simu- 
lations do not describe the real-time dynamics responsible for the generation 
of a finite ^5. On the other hand, the Euclidean instanton transitions (that 
are reproduced in lattice calculations) become suppressed in the deconfmed 
phase, as suggested by the rapid decrease of topological susceptibility above 
T c observed on the lattice [TT] . Because of this, the instantons may no longer 
provide a major contribution to the current susceptibility in the deconfmed 
phase. If so, this would provide a natural explanation for the observation that 
the lattice results are almost independent of the magnetic field at T = 1.12T C 
[37] . This would also agree with the model calculation of Ref. [B] that found 
the CME current becoming insensitive to the magnetic field at high temper- 
ature due to the depletion of instanton effects. One can test this conjecture 
in lattice calculations by evaluating the susceptibility in fixed sectors of topo- 
logical charge (corresponding to a highly excited QCD vacuum configuration 
that arguably resembles the matter produced in a heavy ion collision). 



4-2 Experimental Observable 

In this subsection we will check the consistency of our formula with Ref. [30] • 
For this purpose we shall reiterate some phenomenological discussions in Ref. [30J . 

Using the calculated current and susceptibility we can express the fluctuation 
observables as a function of the volume V, the magnetic field B, the chiral 



24 



chemical potential /15, etc. That is, 



((cos(A0 
c 



Nl 



N± 



fA</> + )» = -((cos(A</>+ + A0_)» 
V 2 N 2 



/ 



2tt 2 



(5f 



The problem is that this expression at first sight may appear different from 
that discussed previously in Ref. [30J. However below we explain that the 
first term stemming from the CME in fact has essentially the same form as 
discussed in Ref. 1301. 



In describing the real data, it is necessary to take account of the screening 
effect in the z direction (denoted previously as y in Ref. [30]). Because /x 5 
emulates the effect of topologically non-trivial domains in a hot medium, it 
would be a reasonable interpretation that our estimate is a valid answer in 
each domain with a finite volume V ~ p 3 . Here, if the topological domain is 
characterized by a sphaleron excitation, p ~ 1/ (a s T) is the typical sphaleron 
size. 

Now the sphaleron can get excited anywhere in the medium and we should sum 
over all excitations. If an excitation occurs deep in a medium, it is unlikely 
that current fluctuations from the corresponding domain can propagate to 
survive outside because of medium screening. On the other hand a sphaleron 
excitation near the surface can easily escape from the medium, and then it 
should contribute to the total current fluctuations if the correlator is for same 
charges. In the case of opposite charges, however, as is clear from our dis- 
cussions in Sec. HJ fluctuations arise from, for example, an upgoing positive 
charge and a downgoing negative charge. Thus, one of two charges must be 
significantly quenched by the medium. In Ref. [30J a phenomenological ansatz 
for such effects was introduced by the following functions; 



g(b/R,\/R) = 



1 



R-b/2 



, , , dx I dz 

2R 2 J-R+b/2 Jz-(x) 



z + {x) 



rR-b/2 rz+(x) 

dx / 

'-R+b/2 Jz-(x) 

for the correlations with same and opposite charges, respectively, where 



1 rR-b/2 rz+ix) 

h(b/R, X/R) = — dx dz z) £_(z, z) 

R z J-R+b/2 Jz^lx) 



exp 



\z±[x) 



z\/X 



(59) 
(60) 

(61) 



with a phenomenological screening length A and 



z+(x) = —z_(x) 



!r 2 -(x- b/2) 2 -R + b/2 < x <0 
. ~ ~ (62) 

IR 2 - (x + b/2) 2 < x < R - b/2 
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is the surface of matter. Here b is the impact parameter, R is the radius of the 
nucleus. 



We postulate that the superposition of the sphaleron domains with V over 
the whole system geometry amounts to the system volume factor with the 
quenching effect taken into account, that is, V is replaced as 




(63) 



= g(b/R,X/R)R 2 Ar]T (for ±± correlations) 
= h(b/R, X/R)R 2 At]t (for ±=p correlations) 

At this point let us think of the /^-integration. The latter term in Eq. 
does not have any /i 5 dependence, so the /i 5 -integration is trivial leading to 
Jd/i 5 >V(/i5) = 1. The former term, in contrast, is proportional to p\. This 
means that p\ after the integration turns to be a parameter (denoted as p 
here) characterizing the dispersion of the p^ distribution. Consequently we 
have the final expressions for the respective cases with the same and opposite 
charges; 



(cos(A0 ± + A0 ± )» = r] + 



N: 



2 



4vr 4 Vf^J ™ 2tt 2 



(64) 

and ((cos(A0 ± + A^))) given by almost the same with V ++ replaced by 
in the above. Here we note that we should identify ^Xg as 



/io oc TArp , (65) 



where T is the sphaleron rate which is proportional to afT 4 at high T and 
thus TAt represents how much topological excitations occur within a time slice 
At. By the dimensional reason p appears. Then we clearly see that the first 
term in Eq. (IM|) has exactly the same structure as the expression discussed in 
Ref. [30J apart from its overall coefficient once it is integrated over the time; 
At — > J dr. In this way we have established a relation between the formalism 
addressed here and the previous work in Ref. [30J. 



To proceed to more concrete analysis we need to specify A, B(t) and N± as 
a function of b (or the centrality). Although those phenomenological analyses 
are important, we will postpone such investigations and discuss them in a 
separate publication which focuses more on the phenomenology of heavy-ion 
collisions. 
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5 Summary 



In this paper we formulated the Chiral Magnetic Effect in terms of the electric- 
current correlation function in a hot QCD medium. We computed the electric- 
current susceptibility in the presence of both the magnetic field B and the 
chiral chemical potential n^. Because of the presence of a preferred direction 
fixed by B, we found that the longitudinal (parallel to B) susceptibility x 
is greater than the transverse (perpendicular to B) one x X - The difference 
arises from only the Landau zero-mode and is given by a UV-finite expres- 
sion; VTN C f <?/|<?/-B|/(27r 2 ). We also gave an intuitive derivation of our result 
based on the anomaly relation. We checked that our result leads to a satisfac- 
tory agreement with the electric-current susceptibility measured in the lattice 
QCD simulation [37] . 

Although x — x 1 ' has an origin in the anomaly, the expression for x ~ X 1 ' 
shows that this difference is independent of /X5 and thus is not sensitive to the 
real-time topological contents of the QCD matter. Therefore we should identify 
it as a background on top of the CME contribution stemming from the square 
of the induced CME current. Since the CME- induced current is proportional 
to B and ^5, the charge- asymmetry fluctuation relevant to the CME has a 
dependence of B 2 and fif which translates into a dispersion parameter of 
the /^-distribution; on the other hand, the non-CME term is proportional 
to B. It would be an interesting question to check if this expected quadratic 
dependence on B can be seen in a lattice simulation restricted to a particular 
topological sector of QCD, that would correspond to a highly excited vacuum 
configuration resembling the matter produced in a heavy ion collision. 

There is an urgent need at the moment in the quantitative phenomenologi- 
cal approaches to CME that on one hand have a firm theoretical ground and 
on the other hand allow a direct comparison to the experimental data from 
RHIC [35]. We view the computation presented here as a necessary step to- 
wards a construction of such an approach. We will report on our approach 
and on the direct quantitative description of RHIC data in the forthcoming 
separate publication. 
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A Appendix: Induced current from the diagrammatic method 



In Refs. [21] and [51] we have discussed several different derivations of the 
induced current in a magnetic field in the presence of non-zero chirality. Here 
we add an alternative derivation using the propagator Eq. (I3T]) . 

From Eq. (131 p we readily obtain for the induced current along the magnetic 
field, 



/ 

x tr 



7 3 P fc (x) (p 



/i 5 7°7 5 



My 



P k (x) 



(A.l) 



Since 7 3 commutes with Pk{x) we need to evaluate the integral over x of Pk{x) 2 
which equals 1 for k > and (1 + isgn(qfB) / j 1/ ~f 2 )/2 for k = 0. Inserting the 
projection operators T±(p) like in Eq. (l4~3l) and using Eq. (f42l it can be seen 
that we need to evaluate the following two traces; 



tr 
itr 



7 1 7 2 7 3 (^57°7 5 + M / )r ± (p) 



-2/x 5 =f 2\p\ 



(A.2) 
(A.3) 



For k > we only need the former trace, which vanishes after integration over 
p z . The only contribution to the current comes from the latter trace when 
k — 0, which is the lowest Landau level. Hence we obtain after summing over 
Matsubara frequencies, 



(J\\) 



VN V ^ f dp " V ^ 5 + 
c y 2tt J 2n ^ ± 2uj pX 



n F {u pX ) 



(A.4) 



where u p x = y(\p z \ + A/is) 2 + Mj for the Landau zero-mode. Noticing that 
(/i 5 + A \p z | ) / cj pA = A duipx/ &p z and integrating over p z we obtain like in Eq. ([151) 
the known result; (^11)^5 = VN c J2f l}B[i 5 / (2tt 2 ), which is independent of My, 
fi q and T. 
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B Appendix: Current-chirality correlation 



Here we will discuss the correlation of the chirality with the current in the 
direction parallel to the magnetic field. We will show that the magnitude of 
this correlation is very similar to the difference between the longitudinal and 
transverse susceptibility. 

By taking the derivative of Eq. (fToT) with respect to /X5 one can easily find 
the correlation between the chirality and the longitudinal current. Such a 
calculation immediately yields, 



(J\\ ^connected = VTN C ^ 7T~Y > (B-l) 



where the total chirality N$ is the volume integral over the zero- component 
of the axial charge, i.e. N$ = J2f /d 3 x(^/7 7 5 ^/). This correlation function 
shows that in a magnetic field the longitudinal current is correlated with the 
chirality, which is the Chiral Magnetic Effect. The magnitude of this correla- 
tion is very similar to the difference between the longitudinal and transverse 
susceptibility as can be inferred from Eq. (|5Tj) . The only alteration comes from 
the charges qf because the chirality is not accompanied by the charge unlike 
the current. However, (Jy iV 5 ) connected depends on the sign of B, while x|i 5 — X^" 5 
does not. 

In the same way as in Sec. 13.51 we could have also arrived at Eq. fIB.ll) using 
the linear response theory. In that case we start from the anomaly equation 
for massless particles which reads 



dx n V 2tt 2 



Applying the same arguments as in Sec. 13.51 and replacing one Jn with N§ in 
the arguments we can recover Eq. fIB.ip correctly. 

The reason why the current-chirality correlator is so similar to the difference 
between the longitudinal and transverse susceptibility is that both quantities 
follow from the axial anomaly. This can be understood explicitly in view of the 
respective starting points of the linear response derivation, namely, Eq. (1521) 
and Eq. flB~2l . 
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